Mixed state properties of superconducting MgB 2 single crystals 
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We report on measurements of the magnetic moment in superconducting MgB2 single crystals. 
We find (M>H° 2 (0) = 3.2 T, po-ffSz (0) = 14.5 T, 7 = 4.6, fi H c {0) = 0.28 T, and k(T c ) = 4.7. The 
standard Ginzburg-Landau and London model relations lead to a consistent data set and indicate 
that MgB2 is a clean limit superconductor of intermediate coupling strength with very pronounced 
anisotropy effects. 

PACS numbers: 74.25.Ha, 74.60.Ec, 74.70.Ad 



The recent discovery of superconductivity in MgB2 
has attracted a lot of attention. Especially the rather 
high transition temperature of nearly 40 K in such a 
simple compound is of interest for applications, but also 
for an analysis of the physical mechanism leading to su- 
perconductivity. Several experiments indicate a phonon 
mediated s - wave BCS mechanism ]2|, [|. Different mod- 
els are proposed to explain the particular properties of 
MgB2 0, H]. Their correctness has to be checked by ex- 
periments, but only a few results are available on single 
crystals g g §, |, 0, P p| . 

We report in this Letter on magnetization measure- 
ments on single crystalline MgB2 in magnetic fields ap- 
plied parallel and perpendicular to the uniaxial crystal- 
lographic (= c) axis. A detailed evaluation allows us 
to obtain the temperature dependence of the most im- 
portant reversible mixed state parameters, such as the 
critical magnetic fields, the characteristic lengths, the 
Ginzburg-Landau (GL) parameter and the anisotropy. 
We will show that MgB2 is a clean limit superconductor 
of intermediate coupling strength with very pronounced 
anisotropy effects. 

Several single crystals of MgB2 were grown using 
high pressure cubic anvils. Details of the process will 
be published elsewhere Two crystals (sample A: 

dx!)xc = 660 x 570 x 21 /j,m 3 ; sample B: a x b x c = 
600 x 384 x 54 /itm 3 ) were investigated by magnetic meth- 
ods. The transition temperature (T c ) of each sample 
was obtained from the ac - susceptibility measured in a 
1 T quantum interference device (SQUID) magnetome- 
ter. Sample A shows an onset of T c at 38 K and a rather 
broad transition of about 1 K. A linear fit of H£ 2 vs. 
T near T c indicates a "bulk transition temperature" of 
37.5 K (see inset of fig. 0a). In sample B we find T c 
= 38.3 K, A T c = 0.3 K and a "bulk T c " of 38.2 K. A 
simple analysis jl4| of the slope of the magnetic moment 
after reversing the applied field demonstrates that the 



size of the domain, in which the supercurrents flow with- 
out impedance, is identical to the sample size. Further- 
more, a comparison of the calculated and the measured 
magnetization in the Meissner regime indicates a super- 
conducting volume fraction of about 100 %. The further 
evaluation of the mixed state parameters did not show 
significant differences between these two crystals. 

The measurements of the magnetic moment were car- 
ried out in the 1 T and in an 8 T (SHE) SQUID mag- 
netometer (for details, cf. (l5)). Fig. [l]a shows the up- 
per critical field of MgB2 for applied fields H a \\c (H^ 2 ) 
and H a \\ab {H^ 2 ). H c2 (T) is determined either from the 
onset of the superconducting signal in the m(T) curve 
("T c (£f a )") or from the disappearance of the supercon- 
ducting signal in the m(H a ) curve. The same results 
were obtained by both methods. H^ 2 could be evaluated 
directly only below 8 T (T > 21 K in this case). At 
lower temperatures the London theory for the reversible 
magnetic moment m r , i.e. m r oc \n{H C 2/H a ), was used 
for the sake of simplicity to extrapolate the experimental 
m r data to zero. The very small magnetic moment in 
higher fields and the logarithmic behavior lead to rather 
large uncertainties in the evaluation, which are indicated 
by error bars in fig. 0a. 

To obtain H^ 2 (0), the data were fitted to the func- 
tion H c2 (t) = H C2 {0){1 - t a f with Moff c c 2 (0) = 3.18 T 
(t = T/T c , T c denotes the bulk transition temperature; 
cv, (3 and H c2 (0) are fit parameters). The initial slope 
of the upper critical field (k = Ho[dH c2 /dT]T c ) is found 
to be -0.112 T/K near T c , thus n H^ 2 (0)/(kT c ) = -0.75. 
This is close to the weak coupling BCS result (= -0.73 
in the clean limit, -0.69 would correspond to the dirty 
limit |l6|]). H c2 (T) is not very sensitive to the coupling 
strength and the above result, -0.75, is actually close 
to that of the strong coupling superconductor Pb [17|, 
of course without considering anisotropy effects |l| 18 1. 
If we apply the same procedure to H a \\ab, we obtain 
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FIG. 1: (a) Upper critical field for H a \\c (# c c 2 ) and H a \\ab 
(H^)- H®2 is obtained from: (i) a direct evaluation of m(T) 
for T > 21 K (solid circles) and (ii) H$ = jH% 2 for T < 21 K 
(solid squares) . The error bars indicate the extrapolation un- 
certainties of the reversible moments measured up to 8 T. The 
BCS curves according to |Ti| are fitted to the experimental 
slope of H c2 near T c . (b) Anisotropy from SQUID {H^/H c c2 ) 
and torque measurements. 



/io_ff"|(0) = 15 T. However, the slope was determined 
in this case from the linear region above 1 T (k = - 
0.55 T/K) and the corresponding extrapolated "T c " of 
36.1 K, because a strong positive curvature of H^iT) is 
observed in the vicinity of T c (cf. the inset of fig. [T]a), 
which represents a well known feature of anisotropic su- 
perconductors, e.g. of the high T' c s, but also of con- 
ventional superconductors, such as Nb. The Fermi sur- 
face and the electron phonon coupling are usually held 
responsible for these anisotropic properties. According 
to jjgfl , the high temperature end of H^(T) can be ob- 
tained by an anisotropic Fermi surface, but not by the 
coupling alone. In Nb, e.g., H C 2(T) is well explained by 
an anisotropic Fermi velocity and an anisotropic electron 
phonon coupling jH| . A similar model could possibly ap- 
ply in the case of MgB2, but alternative theories (e.g. a 
two band model ]^]) are also under discussion. 

The upper critical field anisotropy 7 = H^/H^ is 
shown in Fig. |l|b (full squares). It increases from about 
1 near T c to 4.2 at 22 K, in qualitative agreement with 



previous results 0,0,0. The open squares refer to re- 
sults from torque measurements taken in a 9 T (Quantum 
Design) PPMS system. In this case, the angular depen- 
dence of the reversible torque is fitted to the anisotropic 
London theory with three fit parameters 7, H^ 2 and X a b 
(cf. |pf][). A comparison of the latter two parameters 
with results from the SQUID measurements shows the 
high reliability of the evaluation. Note that this method 
does not lead to the anisotropy of the upper critical field, 
but rather to that of the magnetic penetration depth 
(A), which can, in general, deviate from H^/H^. In 
MgB2, both seem to be the same, at least for T < 30 K. 
The torque indicates a small increase of 7 from 4.3 at 
20 K to about 4.5 at 5 K leading to 7(0) = 4.55. How- 
ever, we cannot exclude some small systematic errors in 
the evaluation, because (i) most of the recorded torque 
data refer to the irreversible regime. Therefore, the re- 
versible signal (t) has to be calculated from the irre- 
versible branches at increasing (t+) and decreasing (r_) 
angles (r = [r+ + r_]/2). The difference between the 
two branches is rather small, but grows at lower temper- 
atures, (ii) The angular dependence of the background 
signal varies with temperature and cannot be determined 
exactly from measurements without a sample. However, 
different data sets for the background do not change 7 
at 5 K significantly. Furthermore, the torque data were 
evaluated at several magnetic fields (0.5 - 2 T), the dif- 
ferences in 7 were very small (2 %). Based on the excel- 
lent agreement between the SQUID and the torque data 
in the overlapping temperature range, we assume that 
A c /A a6 = H^/H c c2 for T < 21 K, which allows us to 
calculate in this temperature range (cf. the solid 
squares for H$ at T < 21 K in fig. |a). This leads to 
Hg$(Q) = 14.5 T. 

The further mixed state parameters can be calculated 
from the London theory and some Ginzburg-Landau re- 
lations. For instance, the magnetic penetration depth 
in the planes (X a b) is obtained from dM/d\n(H a ) = 
</> /(87rA 2 b ) for if a ||c (M = m r /volume, 4> = 2.07 x 
10~ 15 Vs). Since sample A shows a reversible magneti- 
zation already at very small fields (cf. fig. [|) , X a b can be 
evaluated in the whole temperature range from 5 K to T c 
(see fig. ||a). A fit of A~ 6 2 (i) leads to X a b(0) = 82 nm and 
shows that the temperature dependence lies in between 
the (clean limit) BCS Q and a typical strong coupling 
model |jl7). The deviation at lower temperatures indi- 
cates a smaller energy gap - to - T c ratio than according 
to the BCS theory, in agreement with other experiments 
(e.g. |22|, 23 ), which can be explained by the two band 
(assuming a small and a large gap pi|) as well as by the 
anisotropic gap model Further work on more com- 
prehensive calculations including material dependent pa- 
rameters are currently under way. The penetration depth 
in c - direction is obtained from A c = "fX a b, hence A c (0) = 
370 nm. The evaluation of A c from the m(T) measure- 
ments for _ff a ||a6 confirms the anisotropy, but is affected 
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FIG. 2: Temperature dependence of (a) the magnetic pene- 
tration depths, (b) the coherence lengths, (c) the lower critical 
fields, (d) the thermodynamic critical field and the deviation 
function (inset) and (e, f) the GL parameters of sample A. 



by comparatively large errors. 

Further, the (GL) relation hqH"^ = 0o/(27r£„ 6 ) gives 
access to the coherence length in the ab - plane £ a fc and 
in the c - direction £ c = £ Q b/7 (fig- H^)- Accordingly, 
£ ab (0) = 10.2 nm and f c (0) = 2.3 nm. 

The lower critical field H c i can be calculated from 
Mo -ff c x i = [<h,/{2n6\l b )] ■ [HSX ab /U) + 0.5] (a: = c and S 
= 1 for if Q ||c and x = ab and 5 = 7 for H a \\ab), leading 

to /xo-ffci( ) = 63 mT and Mo-tfci( ) = 22 mT (fig. |c). 
A direct experimental assessment of H c \ is usually quite 
difficult, because only the penetration field H Pl i.e. the 
field, at which the first flux lines enter the sample, can 
be obtained, which depends on the sample geometry [ p5[ , 
the anisotropy and the pinning force. We determined H p 
from measurements of the trapped magnetic moment, i.e. 
by measuring the moment in zero field after successively 
applying higher external fields and searching for the first 



deviation from zero at H n 



This procedure is still 



influenced by a finite critical current density ]27] and by 
geometry effects. For example, HqH^ = 4 mT at 5 K 
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FIG. 3: Hysteresis loops of sample A for (a) i/ ]|c and (b) 
H a \\ab. Inset: Irreversibility line of samples A and B for H a \\c 
and H a \\ab. 



(sample A) can be converted into H c i(H p ) for a rectan- 
gular sample geometry p8[ , which leads to H c i/H p = 
17.7 for 7 = 4.5, i.e. Mo# c c i( 5 K ) - 70 mT . but over- 
estimates H^i because the critical current density is not 
taken into account, and is not too far away from the cal- 
culated result of fig. §3 (62 mT at 5 K). 

Furthermore, the thermodynamic critical field is cal- 
culated from the GL relation /j,qH c = 4>o / (VSw ^ab£,ab) 
and found to be 0.28 T at K. Because Af = fJ, H^/2 
(condensation energy), it can also be obtained by inte- 
grating the reversible magnetization M(H a ), i.e. A/ = 
Ho J c2 MdH a . The reversible magnetic moment is ei- 
ther calculated from the irreversible branches of the mag- 
netization in increasing (m+) and decreasing (m_) fields 
in the fully penetrated state, m r = (m+ + m_)/2 (cf. 
fig. ||) or directly measured. The results of the numerical 
integration are shown in fig. |2|d and denoted by H a \\c - di- 
rect and H a \\ ab - direct, respectively. A comparison with 
the GL results indicates that the London model for the 
magnetic penetration depth and the GL relations for H C 2 
and £ represent excellent solutions for MgLV The max- 
imum difference at low temperatures is less than 2 %. 
To check the influence of uncertainties near H p in the 
direct evaluation (geometrical barrier, flux pinning), we 
replace M(H a ) at 5 K by the separately measured Meiss- 
ner slope at < H a < H c i(l — D) and by a simple log- 
arithmic behavior at H c i(l — D) < H a < H c i, i.e. we 
simulate the behavior of an ellipsoidal sample, where the 
" effective demagnetization factor" D is determined from 
M (H a ) = —H a J (1 — 1?) in the Meissner regime. This pro- 
cedure reduces H c at 5 K and brings the above difference 
to almost zero. 

The deviation function, D(t) = [H c (t) / H c (0)]- [1-f 2 ], 
describing the deviation of H c (t) from the parabolic be- 
havior (two fluid model) and indicating the coupling 
strength in a conventional superconductor, is shown in 
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TABLE I: Summary of mixed state parameters for MgB2 



Motf c c 2 (0) 


3.18 T 




14.5 T 


T c 


38 K 


Motf c c i(0) 


63 mT 


HoH2t(0) 


22 mT 




0.28 T 


Ac(0) 


370 nm 


A a6 (0) 


82 nm 


7(0) 


4.6 


€c(o) 


2.3 nm 


U(0) 


10.2 nm 


7(Tc) 


1 


Kf(0) 


8.1 


«f(0) 


37.1 


k(T c ) 


4.7 



the inset of fig. ||d. The maximum of -0.3 - -0.9 % lies 
in between the weak (~ -3.5 %) and the strong coupling 
result (~ +2.5 % for Pb). Although we have to consider 
evaluation errors, the results indicate a clear deviation 
from the weak coupling model, even if we consider the 
anisotropy (cf. J J), and is consistent with other experi- 
ments (e.g. ||, ||). 

The GL parameter k — A/£ is defined at T — T c . 
At lower temperatures the Maki parameters |2{| K\ = 
H c2 /(V2H C ) and k 2 = [0.5 + 0A3/{dM/dH a ) Hc2 - 
0.43D] 1 / 2 can be used with «i(T c ) = k 2 {T c ) = k. ki 
(= A/£ in the GL model) is shown in fig. ^e for H a \\c. Lin- 
ear extrapolations lead to «i(0) = 8.1 and n\ (T c ) = 4.7. 
The ratio k,\{0) / k\{T c ) is 1.72 and considerably larger 
than the BCS value (1.26 in the clean and 1.20 in the 
dirty limit ]l6|| ), but this is not unexpected considering 
stronger coupling JL7| and anisotropy. k 2 depends on the 
slope of M near H c2 and allows a precise determination 
of k c , which is again found to be 4.7 from a linear ex- 
trapolation to T c . For H a \\ab, we get Ki b — jKi and 
therefore nf (0) = 37.1 and «f b (T c ) = 4.7. The errors in 
K 2 b are relatively large in this case because of the very 
small slope of the magnetization near H c2l but the ex- 
trapolation leads to n ab = 5, very close to k c . 

At last, we turn to the irreversible properties of the 
MgB2 single crystals. Hysteresis curves recorded at dif- 
ferent temperatures are presented in fig. |^ for H a \\c and 
H a \\ab. They demonstrate the excellent crystal quality 
by the small hysteresis and the low irreversibility fields 
in both directions. Note that all data points presented 
in fig. [| were measured in the fully penetrated state. 
According to the Bean model J3(| (J c is assumed to 
be constant), the critical current density in the planes 
can be calculated from the irreversible magnetic moment 
(m, = [m+ — m_]/2). For rectangular samples we use 
J C (B) = { mi (B)/n}{4/[b(l - b/3a)}} (sample volume: 
fl = a-b-c), and get 1.4x 10 9 Am -2 at 5 K in the remnant 
state for both samples. To obtain the irreversibility line, 
the onsets of a difference between the field cooled and 
the zero field cooled m(T) measurement were evaluated. 
The results of fig. || (inset) show that the irreversibility 
line is very low for both field directions. 

In summary, we presented measurements of the mag- 
netic moments in single crystalline MgB2 for fields H a \\c 
and H a \\ab, and the subsequent evaluation of the basic 
mixed state parameters. The most important results are 



summarized in table B The general consistency of the 
data set, which is documented nicely, e.g., by the re- 
sults on the thermodynamic critical field, suggests that 
the standard theoretical description can be employed in 
MgB2. The data indicate that MgB2 is a low - k type II 
superconductor in the clean limit with an intermediate 
electron phonon coupling strength (cf . also J3l], [52| ) , but 
a very large anisotropy. 
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